Equilibrium theory of two-dimensional simple liquids

Using the Wigner-Kirkwood expansion and bare Lennard-Jones (LJ) (12-6) potential, an effectiveLJ potential is derived, which includes the quantum effects through the expressions of the effective diameter∂(T, λ) and well-depth (T, λ). We use theWCA perturbation theory to calculate the free energy and pressure for theLJ and effectiveLJ potentials. Simple analytic expressions are given for the reference system and the first order correction calculated. The results are quite good at high density. The quantum effects on the free energy and pressure are also discussed.     

Ab initio quasirelativistic calculations on electronic transitions in ICl by the multireference many‐body perturbation theory

A quasirelativistic perturbative method of ab initio calculations on ground and excited molecular electronic states and transition properties within the relativistic effective core potential approximation is presented and discussed. The method is based on the construction of a state‐selective many‐electron effective Hamiltonian in the model space spanned by an appropriate set of Slater determinants by means of the second‐order many‐body multireference perturbation theory. The neglect of effective spin–orbit interactions outside of the model space allows the exploitation of relatively high nonrelativistic symmetry during the evaluation of perturbative corrections and therefore dramatic reduction of the cost of computations without any contraction of the model‐space functions. One‐electron transition properties are evaluated via the perturbative construction of spin‐free transition density matrices. Illustrative calculations on the X0⁺ ? A1, B0⁺, and (ii)1 transitions in the ICl molecule are reported. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Implementation of the finite field perturbation method in the CRYSTAL program for calculating the dielectric constant of periodic systems

The finite field approach has been implemented in the periodic ab initio CRYSTAL program and been used for calculating the dielectric constants of crystalline LiF and MgO (FCC structure) and BeO (wurtzite structure). To maintain the periodicity along the applied field direction, a "sawtooth" potential is used in conjunction with a supercell scheme. Supercells four to five times longer than the primitive cell in the direction of the applied field provide well-converged results. The influence of the computational parameters is discussed. An alternative scheme has also been implemented, for inner check, that consists of applying a static electric field to a slab of increasing thickness in the direction orthogonal to the surface; the dielectric response at the center of the slab is shown to converge rapidly to the bulk value evaluated with the sawtooth field. The method is accurate and permits the determination of nonlinear corrections to the dielectric constant. When used in conjunction with the local density approximation (LDA) scheme, it provides for the dielectric constant of the three above-mentioned compounds values close to those recently obtained with a time-dependent density functional theory approach.     

Correlated intermolecular interaction components from asymptotically corrected Kohn-Sham orbitals

In these years there was considerable interest inunderstanding of intermolecular forces in energetic(explosive) systems[1—3]. The supermolecular approach(SM) is widely adopted for calculating ab initio in-termolecular interactions. Nevertheless, it is unable toprovide physically meaningful interaction contribu-tions such as electrostatic, induction, repulsion anddispersion energies. In contrast, the symmetry-adaptedperturbation theory (SAPT)[4—8] has the ability to de-rive these correlated…     

SCF perturbation theory for unimolecular reactions

A CNDO/2 SCF perturbation theory is presented for interpreting the form of CNDO/2 potential energy surfaces of unimolecular reactions. The analysis is performed by calculating the energy change E arising from a distortion of the molecular geometry along the reaction coordinate. E is decomposed into different perturbational contributions which are appropriate for an interpretation of the perturbation energy E. Moreover, E is resolved into energy parts arising from a single occupied orbital and contributions due to pairwise orbital interactions. In this way one evaluates numerically how the form of the occupied and unoccupied orbitals determines the magnitude of E. If the distortion occurs along a definite symmetry coordinate, group-theoretical arguments can be applied to discuss the magnitude of characteristic components of the perturbation energy. The SCF perturbation theory is used to analyze the isomerization of ethylene, cis-2-butene and cis-2-butenenitrile.This work was partially supported by Nato-Grant No. 1072     

Improving the efficiency and reliability of free energy perturbation calculations using overlap sampling methods

A challenge in free energy calculation for complex molecular systems by computer simulation is to obtain a reliable estimate within feasible computational time. In this study, we suggest an answer to this challenge by exploring a simple method, overlap sampling (OS), for producing reliable free-energy results in an efficient way. The formalism of the OS method is based on ensuring sampling of important overlapping phase space during perturbation calculations. This technique samples both forward and reverse free energy perturbation (FEP) to improve the free-energy calculation. It considers the asymmetry of the FEP calculation and features an ability to optimize both the precision and the accuracy of the measurement without affecting the simulation process itself. The OS method is tested at two optimization levels: no optimization (simple OS), and full optimization (equivalent to Bennett's method), and compared to conventional FEP techniques, including the widely used direct FEP averaging method, on three alchemical mutation systems: (a) an anion transformation in water solution, (b) mutation between methanol and ethane, and (c) alchemical change of an adenosine molecule. It is consistently shown that the reliability of free-energy estimates can be greatly improved using the OS techniques at both optimization levels, while the performance of Bennett's method is particularly striking. In addition, the efficiency of a calculation can be significantly improved because the method is able to (a) converge to the right answer quickly, and (b) work for large perturbations. The basic two-stage OS method can be extended to admit additional stages, if needed. We suggest that the OS method can be used as a general perturbation technique for computing free energy differences in molecular simulations.     

Direct calculation of ionization potentials of closed-shell atoms and molecules

Vertical ionization potentials, electron affinities and information about quasi-particles can be obtained by using the technique of the single-particle propagator. The expansion of the self-energy part up to third order perturbation theory can be evaluated numerically, but does not lead, in most cases, to satisfying results. A theoretical and numerical analysis of the diagrammatic expansion of the self-energy part requires the introduction of a renormalized interaction and renormalized hole and particle lines.     

Theoretical study of molecular conduction: I. Effective Green's function based on perturbation theory

There have been many experimental and theoretical studies on molecular conduction, as it is a fundamental parameter in the study of molecular‐scale electronics. We have investigated the features of molecular conduction using a Green's function method, which has often been used to solve problems in quantum transport and is also effective in elucidating electron transport in molecules. We have obtained the novel effective Green's functions, including the first‐order energy corrections, by accommodating the self‐energy of the electrodes as perturbation terms. Although these approximate Green's functions only provide information on the first‐order energy corrections, they can involve the elementary properties of molecular conduction. We propose a scheme for the analysis of the relations between molecular orbitals and their roles in molecular conduction and present analytical calculations for normal and cyclic polyenes. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006